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First Results from KamLAND: Evidence for Reactor Anti-Neutrino Disappearance 
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KamLAND has been used to measure the flux of VeS from distant nuclear reactors. In an exposure 
of 162 ton-yr (145.1 days) the ratio of the number of observed inverse /3-decay events to the expected 
number of events without disappearance is 0.611 ± 0.085(stat) ± 0.041(syst) for Ve energies > 3.4 
MeV. The deficit of events is inconsistent with the expected rate for standard Ve propagation at 
the 99.95% confidence level. In the context of two-flavor neutrino oscillations with CPT invariance, 
these results exclude all oscillation solutions but the 'Large Mixing Angle' solution to the solar 
neutrino problem using reactor sources. 

PACS numbers: 14.60.Pq, 26.65.+t, 28.50.Hw, 91.65.Dt 



The primary goal of the Kamioka Liquid scintillator 
Anti-Neutrino Detector (KamLAND) experiment jlj is 
a search for the oscillation of v^s emitted from distant 
power reactors. The long baseline, typically 180 km, en- 
ables KamLAND to address the oscillation solution of the 
'solar neutrino problem' using reactor anti-neutrinos un- 
der laboratory conditions. The inverse /9-decay reaction, 
Ve + p ^ + 11, is utilized to detect D^'s with energies 
above 1.8 MeV in liquid scintillator (LS) [||. The de- 
tection of the e+ and the 2.2 MeV 7-ray from neutron 
capture on a proton in delayed coincidence is a powerful 
tool for reducing background. This letter presents the 
first results from an analysis of 162 ton-yr of the reactor 
j>e data. 

KamLAND is located at the site of the earlier 
Kamiokande |^], with an average rock overburden of 



2,700 m.w.e. resulting in 0.34 Hz of cosmic-ray muons 
in the detector volume. As shown in Fig. |l], the neutrino 
detector/target is 1 kton of ultra-pure LS contained in 
a 13-m-diameter spherical balloon made of 135-/im-thick 
transparent nylon/EVOH (Ethylene vinyl alcohol copoly- 
mer) composite film. The balloon is supported and con- 
strained by a network of kevlar ropes. The LS is 80% 
dodecane, 20% pseudocumene (1,2,4-Trimethylbenzene), 
and 1.52 g/liter of PPG (2,5-Diphenyloxazole) as a fluor. 
A buffer of dodecane and isoparaffin oils between the bal- 
loon and an 18-m-diameter spherical stainless-steel con- 
tainment vessel shields the LS from external radiation. 
During the filling procedure a water extraction and ni- 
trogen bubbling method Q|, optimized for KamLAND, 
was used to purify the LS and buffer oil; PPO prepurifi- 
cation and dust removal were especially important. 
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The specific gravity of the buffer oil is adjusted to be 
0.04% lower than that of the LS. An array of 1,879 pho- 
tomultiplier tubes (PMTs), mounted on the inner sur- 
face of the containment vessel, completes the inner de- 
tector (ID) system. This array includes 1,325 specially 
developed fast PMTs with 17-inch-diameter photocath- 
odes and 554 older Kamiokande 20-inch PMTs . While 
the total photo-cathode coverage is 34%, only the 17- 
inch PMTs corresponding to 22% coverage are used for 
the analysis in this letter. A 3-mm-thick acrylic bar- 
rier at 16.6-m diameter helps prevent radon emanating 
from PMT glass from entering the LS. The containment 
vessel is surrounded by a 3.2 kton water-Cherenkov de- 
tector with 225 20-inch PMTs. This outer detector (OD) 
absorbs 7-rays and neutrons from the surrounding rock 
and provides a tag for cosmic-ray muons. The primary 
ID trigger threshold is set at 200 PMT hits, correspond- 
ing to about 0.7 MeV. This threshold is lowered to 120 
hits for 1 msec after the primary trigger to observe lower 
energy delayed activity. The OD trigger threshold is set 
to provide > 99% tagging efficiency. 

The inner detector is calibrated with 7-ray sources of 
^*Ge, ^^Zn, ^"Co, and Am-Be, deployed at various posi- 
tions along the vertical axis. These sources provide cal- 
ibration energies in the 0.5 to 7.6 MeV region. Energy 
measurements are based on the number of detected pho- 
toelectrons (p.e.) with corrections for PMT gain vari- 
ation, solid angle, the density of PMTs, shadowing by 
suspension ropes, and light transparencies of the LS and 
buffer oil. At the center, about 300 p.e. per MeV are 
observed. Fig. || (a) shows the fractional deviation of 
the reconstructed energies from the source energies. The 
^^Ge and ^°Co sources emit two coincident 7-rays and are 
plotted at an average energy in Fig. ^ (a). The observed 
energy resolution is '-^1 E{Me\f). 

The energy calibration is augmented with studies of 
and ^osrpj (which are contaminants in the detec- 
tor), Bi-Po sequential decays, the spallation products ^^B 



FIG. 2: (a) The fractional difference of the reconstructed av- 
erage 7 energies and average source energies. The dashed line 
shows the systematic error, (b) The B? vertex distribution of 
2.2 MeV neutron capture 7's. 



and ^^N, and 7's from thermal neutron captures on pro- 
tons and ^"^C The reconstructed energy varies by less 
than 0.5% within the 10-m-diameter fiducial volume; lo- 
cal variations near the chimney region are about 1.6%. 
The energy scale exhibits less than 0.6% variation in time 
during the entire data run. Corrections for quenching 
and Cherenkov light production are included, and con- 
tribute to the systematic error in Fig. I (a). The esti- 
mated systematic error for the energy scale is 1.9% at the 
2.6 MeV energy threshold which corresponds to a 2.1% 
uncertainty in the detected neutrino rate. Correlated de- 
cays from Bi-Po and *He/^Li yield energies between 0.6 
MeV and ~15 MeV and are used to extend the range of 
the energy calibration. 

The source positions are reconstructed from the rela- 
tive times of PMT hits. Energy-dependent radial adjust- 
ments are used to reproduce the known source positions 
to ~5 cm; the typical position reconstruction resolution 
is 25 cm. Vertex reconstruction performance through- 
out the LS volume is verified by reproducing the uniform 
distribution of 2.2 MeV capture 7's from spallation neu- 
trons, as shown in Fig. || (b). 

The data presented in this letter were obtained during 
a period from March 4 through October 6, 2002. In to- 
tal we collected 370 million events in 145.1 days of live 
time, corresponding to an average trigger rate of ~30 
Hz. Events with fewer than 10,000 p.e. (-30 MeV) and 
no prompt OD veto signal are classified as 'reactor Vf. 
candidates'; more energetic events are treated as 'muon 
candidates.' 

The criteria for selection of Vf. events are (1) a fidu- 
cial volume cut (i? < 5 m), (2) a time correlation cut 
(0.5 /zsec < AT < 660 /xsec), (3) a vertex correlation 
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cut (Ai? < 1.6 m), (4) a delayed energy window cut (1.8 
MeV < Ed^iay < 2.6 MeV), and (5) a cut on the de- 
layed vertex position requiring it to be more than 1.2 m 
from the central vertical axis (to eliminate background 
from thermometers deployed in this region to monitor 
LS temperature). The overall efficiency associated with 
criteria (2)-(5) combined with the effect of (1) on the 
delayedvertex is (78.3 ± 1.6)%. 

Positrons annihilate in the LS and the detected en- 
ergy includes the kinetic energy plus twice the electron 
rest mass; thus from i>e events produce Eprompt = 
Ejy^ — En — 0.8 MeV, where En is the average neutron 
recoil energy. Anti-neutrinos emitted by ^ssy ^•^^Th 
decays in the Earth, 'geo-neutrinos' [vgeo), contribute 
low-energy events with Eprompt < 2.49 MeV. For exam- 
ple, model la of Dg^o in [p| predicts about 9 Dg^o events 
in our data set. However the abundances of U and Th 
and their distributions in the Earth are not well known. 
To avoid ambiguities from I'geo's we employ (6) a prompt 
energy cut, Eprompt > 2.6 MeV in the present analysis. 

Low-energy 7-rays from ^osq^j j^^g jj^) entering from 
outside the balloon are potential sources of background 
in the range up to 3 MeV and are strongly suppressed by 
the fiducial volume cut (1). The fiducial volume is esti- 
mated using the uniform distribution of spallation neu- 
tron capture events shown in Fig. ^ (b). The ratio of 
these events in the fiducial volume to those in the to- 
tal volume agrees with the geometric fiducial fraction to 
within 4.06%. The same method is used for higher en- 
ergy events, "'^^N, ^^B f3's following muon spallation, and 
agrees within 3.5%. Accounting for uncertainty in the 
LS total mass of 2.1%, we estimate the total systematic 
uncertainty of the fiducial volume to be 4.6%. The den- 
sity of the LS is 0.780 g/cm^ at 11.5°C; the hydrogen- 
to-carbon ratio is computed from the LS components to 
be 1.969 and was verified by elemental analysis (±2%). 
The specific gravity of the LS is measured to 0.01% preci- 
sion using a commercial density meter, and we assign an 
additional 0.1% systematic error due to the uncertainty 
in the LS temperature. The 408 ton fiducial mass thus 
contains 3.46 x 10'^^ free target protons. 

The trigger efRciency was determined to be 99.98% us- 
ing LED light sources. The combined efficiency of the 
electronics, data acquisition, and event reconstruction 
was carefully studied using time distributions of uncor- 
related events from calibration 7 sources. We find that 
this combined efficiency is greater than 99.98%. The ver- 
tex fitter yields > 99.9% efRciency within 2 m of known 
source positions. Using calibrated ^°Co and ^^Zn sources, 
the overall efficiency was verified to within the 3% uncer- 
tainty in source strengths. The detection efficiency of 
delayed events from the Am-Be source (4.4 MeV prompt 
7 and 2.2 MeV delayed neutron capture 7 within 1.6 m) 
was verified with 1% uncertainty. 

From studies of Bi-Po sequential decays, the equi- 
librium concentrations of ^38^ g^j^j ^32^^]^ ^j^g j^g 

are estimated to be (3.5 ± 0.5) x 10~^^ g/g and 
(5.2 ± 0.8) X 10~^^ g/g, respectively. The observed back- 



ground energy spectrum constrains the contamina- 
tion to be less than 2.7 x 10~^^ g/g. The accidental back- 
ground, obtained from the observed flat distribution in 
the delayed time window 0.020-20 sec, is 0.0086 ± 0.0005 
events for the present data set. 

TABLE I: Background summary. 



Background 


Number of events 


Accidental 


0.0086 ± 0.0005 




0.94 ±0.85 


Fast neutron 


< 0.5 


Total B.G. events 


0.95 ±0.99 



At higher energies, the background is dominated by 
spallation products from energetic muons. We observe 
^3,000 neutron events/day/kton. We also expect ~1,300 
events/day /kton Q for various unstable products. 

Single neutrons are efficiently suppressed by a 2-msec 
veto following a muon, but care is required to avoid 
neutrons which mimic the delayed coincidence sig- 
nal. Most fast neutrons are produced by energetic muons 
which pass through both the OD and the surrounding 
rock. This background is evaluated by detecting delayed 
coincidence events with a prompt signal associated with 
a muon detected only by the OD. As expected, a clear 
concentration of events near the balloon edge is observed. 
The expected background inside the fiducial volume is es- 
timated by extrapolating the vertex position distribution 
and considering an OD reconstruction efficiency of 92%. 
To estimate the number of background events due to neu- 
trons from the surrounding rock, the OD-associated rate 
is scaled by the relative neutron production and neutron 
shielding properties of the relevant materials. We esti- 
mate that the total fast neutron background is less than 
0.5 events for the entire data set. 

Most radioactive spallation products simply beta de- 
cay, and are effectively suppressed by requiring a de- 
layed neutron signal. Delayed neutron emitters like ®He 
{Ti/2 — 119 msec) and ^Li (178 msec) are eliminated by 
two time/geometry cuts: (a) a 2-sec veto is applied for 
the entire fiducial volume following a showering muon 
(those with more than 10^ p.e., ^3 GeV, extra energy 
deposition), (b) for the remaining muons, delayed events 
within 2 sec and 3 m from a muon track are rejected. The 
efficiency of these cuts is calculated from the observed 
correlation of spallation neutrons with muon tracks. The 
remaining ^He and ^Li background is estimated to be 



TABLE II: Estimated systematic uncertainties (%). 



Total LS mass 
Fiducial mass ratio 
Energy threshold 
Efficiency of cuts 
Live time 



Total systematic error 



2.1 Reactor power 2.0 

4.1 Fuel composition 1.0 

2.1 Time lag 0.28 

2.1 f spectra §| 2.5 

0.07 Cross section [pH 0.2 



6.4% 
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FIG. 3: Distribution of z/g candidates with fiducial volume 
cut, time, vertex correlation, and spallation cuts applied. The 
prompt energy corresponds to the positron and the delayed 
energy to the captured neutron. The events within the hor- 
izontal lines bracketing the delayed energy of 2.2 MeV are 
due to thermal neutron capture on protons. The events with 
prompt energy below ~0.7 MeV are obtained from the delayed 
trigger. The one event with delayed energy near 4.95 MeV is 
consistent with theexpected 0.54% fraction from '^■^C(n, 7). 



0.94 ± 0.85. The dead time due to the spallation cuts is 
11.4%. This method is checked by exploiting the time 
distribution of the events after a detected niuon to sepa- 
rate the short-lived spallation-produced activities from 
candidates. The uncorrelated event distribution has a 
characteristic time constant of 1/i?^ — 3 sec, where i?^ is 
the incident muon rate. Spallation products have a much 
shorter time constant (~0.2 sec). These methods agree 
to 3% accuracy. As shown in Tabic | the total number of 
expected background events is 0.95 ±0.99, where the fast 
neutron contribution is included in the error estimate. 

Instantaneous thermal power generation, burn-up and 
fuel exchange records for all Japanese commercial power 
reactors are provided by the power companies. The 
time dependence of the thermal power generation data 
is checked by comparison with the independent records 
of electric power generation. The fission rate for each 
fissile element is calculated from these data, resulting in 
a systematic uncertainty in the Ve flux of less than 1%. 
Averaged over the present live-time period, the relative 
fission yields from the various fuel components are ^ss^j 
. 238u : 239py . 24ip^ ^ 0.568 : 0.078 : 0.297 : 0.057. 
The Pe spectrum per fission and its error (2.48%) are 
taken from the literature H] . These neutrino spectra have 
been tested to a few percent accuracy in previous short- 
baseline reactor experiments ^ |^ . The finite /3-decay 
lifetimes of fission products introduce an additional un- 
certainty of 0.28% to the flux; this is estimated from 
the difference of the total De yield associated with shifting 
the run time by one day. The contribution to the flux 
from Korean reactors is estimated to be (2.46 ± 0.25)% 
from the reported electric power generation rates. Other 
reactors around the world give an average (0.70 ± 0.35)% 
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FIG. 4: The ratio of measured to expected flux from reac- 
tor experiments ||l^. The solid dot is the KamLAND point 
plotted at a flux-weighted average distance (the dot size is 
indicative of the spread in reactor distances). The shaded 
region indicates the range of flux predictions corresponding 
to the 95% C.L. LMA region found in a global analysis of 
the solar neutrino data [H. The dotted curve corresponds 
to sin^ 29 = 0.833 and Am^ = 5.5 x lO"'^ eV^ ||l| and is 
representative of recent best-fit LMA predictions while the 
dashed curve shows the case of small mixing angles (or no 
oscillation) . 



contribution, which is estimated by using reactor spec- 
ifications from the International Nuclear Safety Center 
|lO| . The uncertainties for the event rate calculation are 
summarized in Table The errors from reactors outside 
Japan are included in the table under 'Reactor Power'. 

Although the anti- neutrino flux at the location of Kam- 
LAND is due to many nuclear reactors at a variety of 
distances, the De flux is actually dominated by a few re- 
actors at an average distance of ~180 km. More than 
79% of the computed flux arises from 26 reactors within 
the distance range 138-214 km. One reactor at 88 km 
contributes an additional 6.7% to the flux and the other 
reactors are more than 295 km away. This relatively nar- 
row band of distances implies that for some oscillation 
parameters KamLAND can observe a distortion of the 
9e energy spectrum. 

The flux of anti-neutrinos from a reactor a distance 
L from KamLAND is approximately proportional to the 
thermal power flux Ptu/'^T^L'^, where Ptu is the reactor 
thermal power. The integrated total thermal power flux 
during the measurement live time is 254 Joule/cm^. The 
systematic error assigned to the thermal power is con- 
servatively taken as 2% from the regulatory specification 
for safe reactor operation. The corresponding expected 
number of reactor neutrino events (in the absence of neu- 
trino oscillations) in the flducial volume for this data set 
is 86.8 ±5.6. 

The distribution of prompt and delayed energies for 



5 



O 

c 

CD 
> 




2.6 MeV 

analysis threshold 



KamLAND data 
no oscillation 
best-fit oscillation 
sin^2e= 1.0 
Am^=6.9x 10"^eV- 




20 1- reactor neutrinos 

□ geo neutrinos 
|L /y^ □ accidentals 

10 

5 

25 

20 

15 

10 

5 



2 4 6 

Prompt Energy (MeV) 

FIG. 5: Upper panel: Expected reactor z?e energy spectrum 
with contributions of Ugeo (model la of [^) and accidental 
backround. Lower panel: Energy spectrum of the observed 
prompt events (solid circles with error bars), along with the 
expected no oscillation spectrum (upper histogram, with Vgeo 
and accidentals shown) and best fit (lower histogram) includ- 
ing neutrino oscillations. The shaded band indicates the sys- 
tematic error in the best-fit spectrum. The vertical dashed 
line corresponds to the analysis threshold at 2.6 MeV. 



the present sample before the energy cuts are applied 
is shown in Fig. ||. A clear cluster of events from the 
2.2 MeV capture 7's is observed. One event with delayed 
energy around 5 MeV is consistent with a thermal neu- 
tron capture 7 on ^^C. The space-time correlation of the 
prompt and delayed events is in good agreement with ex- 
pectations, and the observed mean neutron capture time 
is 188 ± 23 fisec. After applying the prompt and delayed 
energy cuts, 54 events remain as the final sample. The 
ratio of the number of observed reactor j/g events to that 
expected in the absence of neutrino oscillations is 



N. 



obs 



Nbg 



N, 



expected 



^ 0.611 ± 0.085(stat) ± 0.041(syst). 



The probability that the KamLAND result is consistent 
with the no disappearance hypothesis is less than 0.05%. 
Fig. ^ shows the ratio of measured to expected flux for 
KamLAND as well as previous reactor experiments as a 
function of the average distance from the source. 

The observed prompt energy spectrum is shown in 
Fig. H The expected positron spectrum with no oscilla- 
tions and the best fit with two-flavor neutrino oscillations 
above the 2.6 MeV threshold are shown. A clear deficit 
of events is observed. The measured spectrum is consis- 
tent (93% confidence) with a distorted spectrum shape as 
expected from neutrino oscillations, but a renormalized 



no-oscillation shape is also consistent at 53% confidence. 

The neutrino oscillation parameter region for two- 
neutrino mixing is shown in Fig. ^. The dark shaded 
area is the LMA region at 95% C.L. derived from [ p^ . 
The shaded region outside the solid line is excluded at 
95% C.L. from the rate analysis with Ax^ > 3.84 and 



(0.611 - i?(sin^ 20, Ato2))2 
O.O852 + 0.0412 
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Here, i?(sin^ 29, Am'^) is the expected ratio with the os- 
cillation parameters. 

The spectrum of the final event sample is then ana- 
lyzed with a maximum likelihood method to obtain the 
optimum set of oscillation parameters with the following 
definition: 

X^ato(sin^26', Am2,7VBGi~2,ai~4) 

-2 log Lshapc(sin2 20, Am^, A^bgi~2, ^1-4) 

+Xbg(^BG1~2) + Xdistortion("l~4), 

where Lshape is the likelihood function of the spectrum in- 
cluding deformations from various parameters. iVBGi~2 
are the estimated number of ^Li and *He backgrounds 
and ai^4 are the parameters for the shape deformation 
coming from energy scale, resolution, spectrum and 
fiducial volume. These parameters are varied to mini- 
mize the 3't each pair of [Am2,sin^0] with a bound 
from xIg(^bgi-2) and Xdistortion("i~4). The best fit 
to the KamLAND data in the physical region yields 
sin^ 20 = 1.0 and Am^ = 6.9 x lO'^ eV^ while the global 
minimum occurs slightly outside of the physical region at 



sin 26 = 1.01 with the same Ar 



These numbers can 



be compared to the best fit LMA values of sin 29 = 0.833 
and Am^ = 5.5 x 10"^ eV^ from The 95% C.L. 

allowed regions from the spectrum shape analysis are 
shown in Fig. |6[ The allowed regions displayed for Kam- 
LAND correspond to < < ^ consistent with the solar 
LMA solution, but for KamLAND the allowed regions in 
f < e < f are identical ||l||. 

Another spectral shape analysis is performed with a 
lower prompt energy threshold of 0.9 MeV in order to 
check the stability of the above result and study the 
sensitivity to i>geo- With this threshold, the total back- 
ground is estimated to be 2.91 ± 1.12 events, most of 
which come from accidental and spallation events. The 
systematic error is 6.0%, which is smaller than that for 
the final event sample due to the absence of an energy 
threshold effect. When the maximum likelihood is cal- 
culated, the Pgeo fiuxes from ^ssy g^^^ 232rpj^ ^j,g treated 
as free parameters. The best fit in this analysis yields 
sin^ 26 = 0.91 and Am^ = 6.9 x lO'^eV^. These results 
and the allowed region of the oscillation parameters are 
in good agreement with the results obtained above. The 
numbers of Dgeo events for the best fit are 4 for ^^Sy g^j^j g 
for Th, which corresponds to ^-^40 TW radiogenic heat 
generation according to model la in Q. However, for the 
same model, Vgeo production powers from to 110 TW 
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FIG. 6: Excluded regions of neutrino oscillation parameters 
for the rateanalysis and allowed regions for the combined rate 
and shape analysis from KamLAND at 95% C.L. At the top 
are the 95% C.L. excluded region from CHOOZ and Palo 
Verde |^] experiments, respectively. The 95% C.L. allowed 
region of the 'Large Mixing Angle' (LMA) solution of solar 
neutrino experiments is also shown. The thick dot in- 
dicates the best fit to the KamLAND data in the physical 
region: sin^ 261 = 1.0 and Am^ = 6.9 x 10~^eV^. All regions 
look identical under <-> (tt/2 — 6) except for the LMA region. 



are still allowed at 95% C.L. with the same oscillation 
parameters. 

If three neutrino generations are considered, the D,, sur- 
vival probability depends on two mixing angles 9i2 and 



6*13. In the region close to the best fit KamLAND solution 
the survival probability is, to a very good approximation, 
given by 



P{ve — > i/g) = cos'* 



713 



1 - sin^ 26*12 sin^ 



'12^ 



with Ami2 ^ A?7i^ from the 2-flavor analysis above. 
The CHOOZ experiment |l5[ established an upper limit 
of sin^ 26*13 < 0.15, or co^ > 0.92. The best fit 
KamLAND result would correspond approximately to 
0.86 < sin^ 26*12 < LO. 

In summary KamLAND has used measurements at 
large distances (^^180 km) to demonstrate, for the first 
time, reactor i?e disappearance at a high confidence level 
(99.95%). Since one expects a negligible reduction of 
flux from the SMA, LOW and VAC solar neutrino solu- 
tions, the LMA region is the only remaining oscillation 
solution consistent with the KamLAND result and CPT 
invariance. The allowed LMA region is further reduced 
by these results. Future measurements with greater sta- 
tistical precision and reduced systematic errors will en- 
able KamLAND to provide a high-precision measurement 
of the neutrino oscillation parameters. 
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